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The adsorption of MgO molecules on a Fe(001) surface was studied using density functional theory 
and projector augmented wave methods. The energetically most favored configurations for different 
adsorption sites considered were identified. The most preferable adsorption geometry is when the 
MgO molecules are parallel to the surface, with Mg in the interstitial site and O in on-top of the 
Fe atom. During the adsorption of subsequent MgO molecules in this geometry, a sharp, non- 
oxidized interface is formed between the MgO adlayer and Fe(001) surface. The adsorption of MgO 
perpendicular to the surface, with oxygen incorporated in the topmost Fe layer is less probable, 
but it may lead to the formation of the FeO layer when stabilized with an excess of oxygen atoms. 

Structural, electronic and magnetic properties of both interface types were examined for the MgO 
coverage from 1/9 to 1 monolayer (ML). Electronic and magnetic properties are sensitive to the 
MgO coverage. For lower coverage of MgO, clear hybridization between the Fe 3d and O 2 p states is 
shown. The average magnetic moment of the surface Fe atoms is reduced with coverage, achieving 
2.78 pb for 1 ML of MgO. 

PACS numbers: 73.20.-r, 73.30.+y, 75.70.-i 


I. INTRODUCTION 

Thin-film multilayer systems have attracted consider¬ 
able research interest for many years due to their tech¬ 
nological application in magnetic tunnel junctions, mag¬ 
netic memories (MRAM) or spin valves used in high- 
density magnetic recording devices. In particular, thin- 
film systems containing Fe and MgO layers have been 
investigated by many authors due to their fundamen¬ 
tal properties, such as enhanced magnetic moments^ 
interlayer exchange coupling ^i 3 huge tunneling mag- 
netoresistance^ (TMR), and perpendicular magnetic 
anisotropy— - — 

A small lattice mismatch between MgO (001) and 
Fe(001) surface unit cells, (dM g o ~ V / 2dF e )/dM g o ~ 4%, 
makes favorable conditions for the epitaxial growth of 
both Fe on MgO (Ref. 1771771 ) and MgO on Fe (Ref. [TTUfTTij ) 
with well defined orientation relations: Fe(001) | [100] | 
MgO(001)j[110]. However, the interfacial electronic and 
chemical structure of the Fe/MgO interface becomes 
more complex when considering such effects as Fe-0 bond 
formation and oxidatio n 14 ! 15 This could strongly reduce 
the Fe spin polarization at the interface and thus it makes 
a particular interface structure essential for the TMR ef¬ 
fect. 1 ® Experimentally, many different, apparently con¬ 
tradicting structures at the Fe/MgO and MgO/Fe in¬ 
terfaces were reported, ranging from nearly iron-oxide- 
freeii - — to the formation of an FeO layer— i - — The ob¬ 
served differences originate from preparation conditions, 
either oxygen-deficient (deposition from bulk MgO) or 
oxygen-rich (reactive deposition of Mg in oxygen atmo¬ 
sphere), as well as the deposition sequence (Fe on the 
MgO surface or MgO monolayers on Fe). In particular, 


the interfaces between the MgO and Fe layers (Fe on MgO 
and MgO on Fe) were experimentally examined by con¬ 
version electron Mossbauer spectroscopy (CEMS)— 4 Un¬ 
der oxygen-deficient conditions, the coexistence of both 
oxidized and nonoxidized interfaces was confirmed, and 
interface properties were found to be highly sensitive to 
the concentration of defects—! Different methods have 
been applied to prevent the formation of FeO layers under 
oxygen-rich conditions, including modified reactive depo¬ 
sition^ deposition of metallic Mg atoms on the Fe(001) 
substrate and subsequent annealing of the sample under 
O 2 exposure^ or a metallic Mg buffer-layer formation 
followed by reactive Mg deposition. 27 On the theory side, 
a full-potential linearized augmented-plane-wave study of 
monolayers of Fe on MgO(OOl) has shown!! that there is 
almost no interaction between Fe and MgO. More recent 
plane-wave DFT calculation s 29 ! 30 of MgO monolayers on 
a Fe(001) substrate confirmed that the MgO films of one 
to three monolayers (MLs) interact only weakly with the 
Fe substrate. At the oxidized MgO/FeO/Fe(OOI) inter¬ 
faces, an increase of the magnetic moments of Fe atoms 
near the interface and a reduction of the work func¬ 
tion with respect to the clean Fe(001) surface were re¬ 
port ed— 

The first stage of growth, i.e. adsorption, plays a signif¬ 
icant role in the structure of multilayer systems. Under¬ 
standing of this initial process allows for an understand¬ 
ing of the nucleation and growth of three-dimensional 
structures. The MgO substrate is commonly used due 
to its simple atomic structure and its relevance to many 
experimental works on model systems— ! Experimental!! 
and theoretical studie s 33 ! 34 of Fe adsorption on defect- 
free MgO(OOl) surface showed the strongest adsorption 
binding for Fe adatoms on-top of the surface oxygen 
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atoms. Surprisingly, there are no theoretical works con¬ 
sidering the adsorption of single MgO molecules on metal 
surfaces. The equilibrium structures of the MgO layers 
on Fe(001) were optimized with an a priori assumption 
that the oxygen preferably adsorbs at the on-top sites. 
However, the calculation by Beltran et suggested 
that under low MgO coverage, Fe binds to Mg rather 
than to O. 

Simulations of MgO adsorption on Fe(001) should ac¬ 
count for different deposition methods. The reactive de¬ 
position of metallic Mg in a molecular oxygen atmosphere 
should be avoided because it leads to direct oxidation 
of the Fe surface layer^ For the preferred deposition 
of MgO from the bulk, one should consider the adsorp¬ 
tion of MgO molecules and/or of the products of their 
dissociations, i.e. co-adsorption of atomic (ionic) oxy¬ 
gen and magnesium. In the former case, one can refer 
to the adsorption of atomic oxygen on iron surfaces^ 
which confirmed the experimental observation 37 that the 
interstitial (fourfold hollow) position is the most stable 
oxygen adsorption site on the Fe(001) surface. Vassent 
et air 8 experimentally demonstrated that when MgO is 
deposited from a bulk target by electron beam bombard¬ 
ment, the resulting beam is composed mainly of atomic 
Mg and O, which recombine on the Fe(001) substrate. 
However, adsorption of the MgO molecules should also 
be considered, especially when an alternative deposition 
method is used (e.g. pulsed laser deposition or thermal 
evaporation). A simple mass-spectroscopy experiment 
performed for the purpose of the present study demon¬ 
strated that when MgO is evaporated from a thermal 
source, the vapor flux, contains, aside from atomic O and 
Mg, at least 15% of the MgO molecules^ The real contri¬ 
bution of the MgO molecules can be considerably higher 
when considering the dissociation of the MgO molecules 
in the quadrupole mass spectrometer. 

In this paper, we report DFT calculations of the ad¬ 
sorption of MgO molecules on the Fe(001) surface. Start¬ 
ing from a single MgO molecule in a large surface unit 
cell and by considering different adsorption sites and 
molecule orientations, the most energetically favorable 
adsorption geometries were determined. Subsequently, 
by increasing the coverage to the complete MgO mono- 
layer, we were able to simulate low-temperature de¬ 
position of MgO onto the clean Fe(001) surface under 
ultrahigh-vacuum conditions. We have examined the en¬ 
ergetics and electronic properties of the ensuing struc¬ 
tures, which are related to the first stages of MgO growth 
on Fe(001). 


II. CALCULATION METHOD 

The calculations presented in this work were car¬ 
ried out using the Vienna ab initio simulation 
package (VASP) based on density functional theory 
(DFT) . 40 ! 41 The exchange-correlation energy was de¬ 
scribed within the Perdew-Burke-Ernzerhof (PBE) ver¬ 


sion 4 -. of the spin-polarized generalized gradient approx¬ 
imation (GGA)i 43 ! 44 The electron-ion interactions were 
represented by the projector augmented-wave (PAW) po¬ 
tentials^ with the 3d 7 4s 1 , 3s 2 , and 2s 2 2p 4 states consid¬ 
ered as the valence states of the Fe, Mg, and O atoms, 
respectively. For the plane-wave-basis set, a cutoff energy 
of 500 eV was applied based on the series of tests. The 
fc-space integrations for the bulk bcc Fe were performed 
using al2xl2xl2 special Appoint mesh generated by 
the Monkhorst-Pack method^ The partial occupancies 
were treated by using the first-order Methfessel-Paxton 
method^ with a Fermi surface smearing of 0.2 eV. This 
computational setup yielded calculated properties of bulk 
bcc iron that were in good agreement with experiment 
and previous theoretical works. The optimized lattice 
constant of 2.832 A for ferromagnetic bcc Fe differs from 
the experimental value (2.866 A)^& by only approxi¬ 
mately 1% and agrees well with the values reported in 
other theoretical works i 30 ! 49 — The bulk modulus, ob¬ 
tained by fitting the Murnaghan equation of state to 
the energy dependence on the lattice parameter, is 179 
GPa, which is slightly larger than the measured value 
of 172 GPa^i and it is close to the other theoretical 
GGA^r— and the local density approximation (LDA)^ 
results. The calculated cohesive energy of 4.943 eV is 
15% larger than the experimental value (4.28 eV)^ but 
it is similar to previous results obtained by other authors 
using both pseudopotential^ and full-potentia l 56 ! 57 first- 
principles methods. The computed magnetic moment of 
2.19 pb is also in good agreement with the experimental 
bulk value of 2.22 and other calculations 4 ^ 53 i 58 

The Fe(001) surface was modeled by a slab of 9 Fe lay¬ 
ers separated from its periodic replicas by a thick vacuum 
region of 22 A. An adequate number of Monkhorst-Pack 
fc-points has been selected for each system depending on 
the size of the unit cell: 12 x 12 x 1 for the clean Fe(001) 
surface lxl cell and 4 x 4 x 1 for MgO adsorption in the 
3x3 surface unit cell. MgO molecules were adsorbed on 
both sides of the slab. The positions of all atoms were 
optimized until the forces exerted on each atom were less 
than 0.01 eV /A. 

The adsorption energy was calculated from the total 
energy difference: 


E&d — Q/Fe(001) — ^Fe(001) — AbE Mg o), (1) 


where E , Mg0/ / Fe ( 001 ) and ^Fe(ooi) represent the energy of 
the slab with adsorbed MgO molecules and the total en¬ 
ergy of the clean slab, respectively. N is the number of 
molecules adsorbed on both sides of the slab, and Em s o 
is the energy of the free MgO molecule calculated in a 
large rectangular box with dimensions of 11 x 12 x 13 
A. The calculated bond length of a free MgO molecule is 
1.75 A. 
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FIG. 1: (color online) Adsorption sites (crosses) for MgO 
molecules in the 3x3 surface unit cell representing the Fe(001) 
surface: ot - on-top, int - interstitial, bri - bridge site. The 
brown and light brown balls correspond to the first and second 
Fe layers, respectively. 

III. RESULTS AND DISCUSSION 
A. Clean Fe(001) surface 

The calculated basic properties of the clean Fe(001) 
surface, lattice relaxation, surface energy, work function, 
and magnetic moments on the surface atoms create a 
solid reference system for MgO adsorption. The ob¬ 
tained values of relaxation, A ij = (dij — d)/d, of surface 
atomic layers distance, dij, between subsequent layers i 
and j, where d is the interplanar distance in bulk Fe, 
fall well (Table HT1) within the range of values determined 
by calculations applying the sam e 30 ' 58 or a similar fla¬ 
vor of GGA-P W9 1, 51 1 59 The computed surface energy, 
E sur i = 2.52 J/m 2 is between two experimental results, 
2.41 J/m 2 (Ref. ]60) and 2.55 J/m 2 (Ref. l6ll ). and is of 
a similar magnitude as the result reported by other au- 
thors i 30 : 51 : 52 This energy was determined from the ex¬ 
pression E SUI { = jj(E n — tiEb), where E n is the total 
energy of the slab, n is the number of layers in the slab, 
A is the surface area, and the 1/2 factor accounts for the 
two surfaces of the slab. Eb is the energy of the bulk 
layer calculated as the difference between the total en¬ 
ergy of the (n +1)- and n-layer slabs. The work function 
of 3.86 eV, calculated as the difference between the elec¬ 
trostatic potential energy in the vacuum region and the 
Fermi energy of the slab, is in perfect agreement with ear¬ 
lier GGA calculations ; 30 ! 51 : 58 However, it is substantially 
less than the experimental values of 4.67 eV (Ref. l6a ) 
and 4.24 eV (Ref. |63). The magnetic moments on sur¬ 
face and subsurface Fe atoms of the Fe(001) slab, due to 
their lower coordination, are enhanced to 2.94 and 2.33 
Hb, respectively. They agree well with previou s 30 ! 51 i 52 i 58 
GGA results. For the deeper layer atoms, the magnetic 
moments approach the value for bulk iron of 2.19 ^b- 

B. MgO molecules on the clean Fe(001) surface 

1. Geometry and energetics 

MgO molecules were adsorbed on the relaxed Fe(001)- 
oriented slab. Figure [T] shows the adsorption sites con- 


TABLE I: Adsorption energy S a d for a single MgO molecule 
adsorbed on the Fe(001) surface for all initially considered 
geometries. For the parallel geometry, the first part of each 
term describes the location of the Mg and the second part 
describes the location of the O atom. For the perpendicular 
geometry, the first part of the term describes the atom that 
bonds with the surface and the second part determines the 
location of the molecule. 


Parallel geometry 

Perpendicular geometry 

Configuration 

S ad [eVj 

Configuration 

S ad [eVj 

int-ot 

4.60 

O-ot 

3.31 

ot-int 

4.37 

O-int 

4.19 

bri-ot 

4.17 

O-bri 

3.83 

ot-bri 

4.14 

Mg-ot 

1.38 

bri-int 

4.40 

Mg-int 

0.64 

int-bri 

4.90 

Mg-bri 

1.33 


sidered for Mg and O atoms of the MgO molecule at the 
Fe(001) surface: on-top (ot) of the Fe atom, the inter¬ 
stitial (int) four-fold hollow between the four surface Fe 
atoms and the bridge (bri) site between the two near¬ 
est surface Fe atoms. Adsorption geometries tested in 
the calculations included different possible alignments of 
the MgO molecules with respect to the surface: parallel, 
with Mg and O at different positions, and perpendicu¬ 
lar, with the Mg or O atom bonding with the substrate. 
This resulted in six orientations of MgO molecules par¬ 
allel to the surface, termed int-ot, ot-int, bri-ot, ot-bri, 
bri-int, and int-bri, where the first part of each term in 
this notation describes the location of the Mg and the 
second part describes the location of the O atom. Sim¬ 
ilarly, six configurations of the MgO molecule oriented 
perpendicular to the surface were considered: three of 
them bonding with O (O-ot, O-int, O-bri) and the re¬ 
maining three bonding with the Mg atom to the Fe(001) 
surface (Mg-ot, Mg-int, Mg-bri). 

The adsorption energy for a single MgO molecule ad¬ 
sorbed on the 3x3 surface unit cell (a coverage of 
1/9 ML) for different adsorption geometries is shown in 
Table Q] Configurations aligned parallel to the surface 
have the highest adsorption energy, especially the int-ot 
and the int-bri systems with values of 4.60 eV and 4.90 
eV, respectively. Generally, the bonding of the single 
MgO molecule adsorbed perpendicularly is substantially 
weaker, especially via the Mg (S a d «1 eV compared to 
less than 4 eV for the binding via the O atom). The only 
exception is the O-int configuration, whose adsorption 
energy of 4.19 eV is comparable with the parallel MgO 
configuration. In this case, the oxygen atom is at the 
closest distance to the surface, 0.69 A, and its local con¬ 
figuration is similar to that in the FeO(OOl) monolayer. 

A coverage dependence of the stability of different ad¬ 
sorbate structures is discussed below. Starting from a 
single MgO molecule on the 3x3 surface unit cell, by the 
subsequent addition of an extra MgO molecule, in all pos¬ 
sible arrangements with respect to the previous coverage 
and preserving the appropriate adsorption configuration, 
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FIG. 2: (color online) Coverage dependence of adsorption en¬ 
ergy (a), Mg-O distance (b) and work function (c) for the 
most favorable overall int-ot orientation of MgO molecules 
adsorbed on the Fe(001) surface. For comparison, the adsorp¬ 
tion energy plot in (a) also presents the results for the most 
favorable perpendicular O-int orientation of the molecules. 



the coverage dependence was followed up to a complete 
MgO monolayer. In this manner, the lowest-energy ar¬ 
rangement for each degree of coverage was determined. 
Only the int-ot configuration with MgO aligned paral¬ 
lel to the surface and the perpendicular O-int configura¬ 
tion led to the experimentally observed structure, i.e. to 
the pseudomorphic Mg0(001)/Fe(001) interface. Other 
parallel configurations of a single MgO molecule, includ¬ 
ing int-bri, which shows the strongest adsorption bind¬ 
ing for 0=1/9 ML, are at higher coverages considerably 
more weakly bound than the int-ot phase, and they do 
not form any reasonable (known) structure. In turn, at 
higher coverages, the binding of MgO in the perpendic¬ 
ular O-int configuration is distinctly weaker than in the 
int-ot. Therefore, and taking into account the highest 
adsorption energy for higher degrees of coverage, we con- 


FIG. 3: (color online) Side and top (right column) views of 
the most stable final int-ot configuration of MgO molecules 
adsorbed on the Fe(OOl) surface, with the Mg atom (red) in 
the interstitial (int) and the O atom (blue) on-top (ot) of the 
Fe atom. Medium (brown) balls represent Fe atoms. In the 
top view, Fe atoms of the second layer are a lighter shade. 


sider the int-ot system as the most stable configuration. 
Therefore, in the following, we concentrate on the paral¬ 
lel int-ot system and some aspects of the perpendicular 
O-int configuration in which the surface Fe layer is prone 
to oxidation. 

The adsorption energy as a function of MgO coverage 
is plotted in Fig. Ha). In the int-ot configuration, with 
the Mg atoms in the interstitial sites and the O atoms in 
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the on-top sites, the adsorption energy increases mono- 
tonically with the coverage, from 4.60 eV for 1/9 ML to 
6.78 eV for 1 ML of MgO. For comparison, in Fig.j^a) we 
have also plotted the variation of the adsorption energy 
versus coverage for MgO molecules adsorbed in the per¬ 
pendicular O-int configuration, with the O atom closer 
to the surface. In this case, the MgO binding is weaker 
than for parallel-oriented MgO molecules. It increases 
only slightly with increasing coverage, and the difference 
in the adsorption bond strength for 1 ML of coverage is 
larger than 2 eV. This shows that when stoichiometric 
MgO is adsorbed, the formation of the interfacial FeO 
layer is less favorable. 

The optimized configurations for several coverages of 
the int-ot-oriented MgO molecule are shown in Fig. [3] As 
can be seen, MgO molecules adsorbed on the surface have 
a tendency to form compact, symmetric two-dimensional 
MgO surface clusters. This is in agreement with the ex¬ 
perimentally observed tendency of island growth during 
MgO deposition^ Any asymmetry in the cluster shape 
would lead to atom displacements from an ideal adsorp¬ 
tion position, as observed for 0=5/9, where one of the O 
atoms is clearly shifted with respect to the other atoms 
due to a smaller number of the nearest neighbours. Fig¬ 
ure [^b) shows the change of the Mg-0 bonding distance 
with the degree of coverage. Upon interaction with the 
surface, (Aig-O increases from the initial 1.75 A in a free 
molecule to 2.01 A for the full ML coverage. This value is 
lower than the value in bulk MgO (2.11 A) ^ which is due 
to the pseudomorphicity with the Fe(001) surface. In the 
int-ot configuration, a minimum for 5/9 ML coverage is 
noted, which results from the previously mentioned lower 
symmetry of the MgO adsorbate structure. The 4% lat¬ 
tice misfit between MgO and the Fe(001) surface causes 
the Mg atoms to be closer to the surface than O atoms by 
0.11 A. Such a rumpling of the MgO layer closely repro¬ 
duces the results of Jeon et al£2 Whereas the complete 
monolayer is pseudomorphic, for incomplete MgO cov¬ 
erage the interaction between Mg and O atoms on the 
surface causes the O atoms not to be located directly on- 
top of Fe atoms but to be slightly shifted toward the Mg 
atoms (Fig. 02). 

The calculated height of the O atom above the Fe sur¬ 
face, ho, increases with the degree of coverage from the 
initial value of 1.90 A for O =1/9 to 2.24 A for 0 =1. Our 
results also agree with the measured Fe-0 distance, 2.20 
± 0.05 A, at the Fe(001)/Mg0(001) interfaced At 1 ML 
coverage, for the int-ot adsorption geometry, there is no 
incorporation of oxygen into the surface Fe layer. This 
configuration represents the sharp nonoxidized MgO/Fe 
interface, as discussed in the so-called Mg-rich condi¬ 
tion^ 

Table [TT] demonstrates that MgO adsorption in the int- 
ot configuration does not cause significant changes in re¬ 
laxation of the interlayer distances at the Fe(001) surface. 
A weak derelaxation can be seen at fractional coverages. 
However, for 1 ML coverage, relaxations are nearly iden¬ 
tical to the clean Fe(001) surface, indicating that the 


TABLE II: Relaxations A ij (in % of the bulk interplanar 
spacing) of the three top interlayer distances for the Fe(001) 
surface as a function of the MgO coverage in the int-ot con¬ 
figuration. Positive and negative values correspond to the 
expansion and contraction of the spacing, respectively. 


Coverage (ML) 

Al2 

A 23 

A 34 

0 

- 2.1 

3.4 

0.7 

1/9 

- 1.2 

2.7 

0.4 

2/9 

- 0.8 

2.6 

0.4 

3/9 

-0.3 

2.1 

0.4 

5/9 

0 

2.1 

0.6 

7/9 

- 0.6 

2.1 

0.6 

9/9 

- 2.1 

2.6 

0.9 


Fe(001) surface is stable against the adsorption of the 
MgO monolayer, when a sharp metal-oxide interface is 
formed. 

The interface structure is highly sensitive to the ini¬ 
tial orientation of the MgO molecule with respect to the 
Fe surface. When the perpendicularly orientated MgO 
molecules approach the surface to form the O-int con¬ 
figuration, the O atom oriented toward the surface is 
incorporated into the Fe surface. For a complete MgO 
monolayer, the O atoms are 0.725 A above the Fe(001) 
surface. In this case, formation of the Fe-0 layer is possi¬ 
ble. This would result in the oxidized MgO/Fe interface 
with an uncompensated surface layer of MgO. To make 
this case more realistic, we simulated O-rich conditions 
by providing additional oxygen atoms adsorbed in the 
on-top positions. This resulted in a stable surface MgO 
monolayer and a subsurface FeO layer. The FeO layer 
showed only a small rumpling of 0.388 A, which corre¬ 
sponds to a nearly perfect oxidized MgO/FeO/Fe(001) 
interface, in agreement with other calculations : 29 ' 67 

The variations in the work function with changes in the 
coverage of MgO adsorbed in the int-ot configuration are 
plotted in Fig. [He). The work function decreases grad¬ 
ually from 3.85 to 2.45 eV. This large reduction results 
from a difference between a small increase in the vac¬ 
uum electrostatic potential and a pronounced shift of the 
Fermi energy of the MgO/Fe(001) system. The changes 
in work function with coverage reflect the changes in the 
electron charge-density distribution induced by MgO ad¬ 
sorption. Figure^] diplays the isosurfaces of the electron 
density difference: 

Ap(r) = p MgO/Fe(001)( r ) _ p Fe(001)( r ) _ p MgO( r ) > 

where ^MsO/FU 001 ) j s the electron density of the Fe(001) 
surface with adsorbed MgO, pFe(ooi) is the electron den¬ 
sity of the Fe(001) surface with the frozen atomic posi¬ 
tions, and p M s° is the electron density of the frozen MgO 
adsorbate. For a single MgO molecule (O = 1/9), the iso¬ 
surfaces show the local electron charge accumulation in 
the space between the Mg and the neighboring Fe atoms. 
This means that MgO binding to the surface is domi¬ 
nated by the Mg-Fe bonding, which favors a horizontal 
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FIG. 4: (color online) Top and side (bottom panels) views of isosurfaces of the valence charge density difference, Ap(r), 
depending on the MgO coverage. Electron charge gain/loss is drawn in yellow/blue. Isosurface density level is 0.005 e/bohr 3 . 
Red, blue, and brown balls represent Mg, O and Fe atoms, respectively. The black lines in the top view panels mark the plane 
cuts for the side views shown in the lower row panels. 
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FIG. 5: (color online) Bader charge difference (A Q) on the 
surface atoms of the top layers of the MgO/Fe(001) system 
induced by different degrees of MgO coverage. 


alignment of the MgO molecule. In general, the plots 
for different degrees of coverage show that the electron 
charge is drawn from the less negative Mg atoms (Paul¬ 
ing electronegativity =1.31) toward the more electroneg¬ 
ative Fe atoms (electronegativity =1.83). This results 
in a negative contribution to the surface dipole barrier 
of the clean Fe(001) and the lowering of the work func¬ 
tion. However, this effect is partially compensated by 
electrons drawn by the strongly electronegative oxygen 
ions (electronegativity =3.44). Consequently, the elec¬ 
trostatic dipole barrier changes relatively little (by ~ 0.4 
eV) in the range of coverage considered, and the work 
function lowering is mainly due to the shift of the Fermi 


level to higher energies with MgO coverage. 

More quantitative information about the charge trans¬ 
fer for each atom is provided by the analysis of the cal¬ 
culated Bader charge^ difference: 

A Q = QMgO/Fe(001) _ QFe(OOl) _ gMgO 

where Q M g°/ Fe (ooi) i s the Bader charge of the atoms of 
the MgO/Fe(001) system, Q Fe ( 001 ) is the charge on the 
atoms of the frozen Fe(001) surface, and Q Mg0 is the 
charge of the frozen MgO molecule. The largest charge 
difference occurs at the surface Fe layer atoms (Fig. [5]). 
The O atoms and Fe atoms located under the oxygen lose 
electron charge, which is transferred to the Mg atoms and 
Fe atoms away from the adsorbate. An exception occurs 
when the O atom has a lower number of neighboring 
atoms in plane (e.g., for 0=1/9 and 0=5/9). In this case, 
oxygen gains electron charge. Much less electron transfer 
can be observed when the MgO monolayer is completed. 
The analysis of the total Bader charge shows transfer 
from oxygen to the Fe surface, which increases with the 
coverage, except for 4/9 and 5/9 ML. It can reduce the 
surface dipole barrier and consequently decrease the work 
function. 


2. Electronic and magnetic structure 

The partial, layer-resolved densities of states (PDOSs) 
for the considered system are presented in Fig. [Gj For the 
clean Fe(001) surface, the changes in the PDOSs with re¬ 
spect to the bulk are practically limited to the atoms of 
the topmost Fe layer. For the first surface Fe layer, the 
PDOSs were plotted separately for the Fe atoms away 
from the O atoms and for those binding with the O 
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FIG. 6: (color online) Partial densities of states (PDOSs) of 
the clean Fe(001) surface and for different degrees of MgO cov¬ 
erage adsorbed in the most stable int-ot configuration. Pos¬ 
itive and negative densities correspond to the spin-up and 
spin-down states, respectively. Energy zero corresponds to 
the Fermi level. The plots for Mg states are multiplied by 5. 
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FIG. 7: Coverage dependence of the average magnetic mo¬ 
ment on the Fe atoms in the topmost surface layer. 
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face induces changes in the magnetic moment of the Fe 
atoms of the topmost surface layer. The values of the lo¬ 
cal magnetic moment vary between 3.00 p,s for Fe atoms 
without contact with an adsorbate and 2.59/is for Fe 
atoms adjacent to the MgO molecules. The average mag¬ 
netic moment of the surface Fe atoms plotted in Fig.[7]as 
a function of MgO coverage gradually decreases with the 
coverage, from 2.94 ps for the clean surface to 2.78 pb for 
1 ML of MgO. The latter value is 5% larger than the value 
previously reported for 1 ML of MgO on Fe(OOl)^ Gen¬ 
erally, the proximity of oxygen enhances the magnetic 
moment of iron^FS 7 which means that nonmagnetic Mg 
atoms are responsible for weakening the magnetism on 
the Fe surface. As a result of the Mg-atom proximity, the 
Fe 3d majority states are shifted toward higher energy. 
Their occupancy decreases, resulting in the reduction of 
the magnetic moment of Fe. The adsorption of MgO has 
a negligible effect on the magnetic moments of the sec¬ 
ond Fe layer, which are nearly unaltered compared to the 
corresponding moments of the clean Fe(001) surface. Hy¬ 
bridization of the Fe 3d states and O 2 p states induces a 
small magnetic moment on the O atoms, which changes 
with coverage. The magnetic moment ranges between 
0.17 pb and 0.03 ps for 0 = 1/9 and 1 ML of coverage, 
respectively. A small magnetic moment of 0.01-0.03/is 
also appears at the Mg atoms. 


atoms. In the entire range of coverage, the former are 
nearly unaltered compared with the Fe 3d states of the 
topmost layer of the clean Fe surface. The main changes 
are observed in the PDOSs of the Fe atoms located under 
the oxygen. This effect is particularly visible for low cov¬ 
erages. The proximity of the O atoms contributes to the 
additional peaks in the Fe 3d states due to hybridization 
with the O 2 p states. The additional peaks disappear 
with the increasing coverage. For the full MgO mono- 
layer, the 3d PDOS of the topmost Fe atoms is similar to 
that for the clean surface atoms, but the PDOS peaks are 
sharper and higher. With increasing coverage, the shape 
of the O 2 p PDOS changes and the states are shifted to 
lower energies, which reflects the stronger MgO binding 
to the iron substrate. 

The presence of the Mg and O adatoms on the sur- 


IV. SUMMARY AND CONCLUSION 

We have presented the results of a DFT study of the 
adsorption of MgO molecules on a Fe(001) surface for 
coverage varying from 1/9 to 1 ML. The most stable 
configuration of the adsorbed MgO molecule is with the 
Mg atoms in the interstitial site and the O atoms on-top 
of Fe site (int-ot configuration) for the molecule aligned 
parallel to the surface. The configuration with Mg at 
the interstitial and O at the bridge site (int-bri) that is 
preferred at the initial adsorption stages becomes unsta¬ 
ble for higher degrees of coverage. This confirms the a 
priori assumption that the int-ot configuration of the ad¬ 
sorbed MgO on the Fe(001) surfac o 26 ’ 69 is the preferred 
arrangement. Under circumstances that simulate the ad- 




























sorption of stoichiometric MgO, we show a preference for 
the sharp MgO/Fe interface formation without oxidation 
of the topmost Fe layer. 

Our results also show that oxidation of the Fe(001) sur¬ 
face cannot be excluded during the adsorption of MgO 
molecules aligned perpendicular to the Fe surface with 
the O atoms facing the substrate. The analysis of the 
adsorption energy shows that this configuration is less 
probable in stoichiometric conditions than the adsorp¬ 
tion of MgO molecules aligned parallel to the surface. 
Nevertheless, in this particular O-int configuration, the 
O atoms significantly approach the surface and an FeO- 
like layer is formed. This is consistent with many exper¬ 
imental works22r—, reporting the appearance of the FeO 
phase at the interface in MgO/Fe systems. 
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